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Abstract. When stars like our Sun are young they rotate rapidly and are very magnetically 
active. We explore dynamo action in rapidly rotating suns with the 3-D MHD anelastic spherical 
harmonic (ASH) code. The magnetic fields built in these dynamos are organized on global-scales 
into wreath-like structures that span the convection zone. Wreath-building dynamos can undergo 
quasi-cyclic reversals of polarity and such behavior is common in the parameter space we have 
been able to explore. These dynamos do not appear to require tachoclines to achieve their spatial 
or temporal organization. Wreath-building dynamos are present to some degree at all rotation 
rates, but are most evident in the more rapidly rotating simulations. 
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1. Introduction 



When stars like the Sun are young, they rotate quite rapidly. Observations of these 
young Suns indicate that they have strong surface magnetic activity and can undergo 
global-scale polarity reversals similar to the 22-year solar cycle. The magnetic fields 
observed at the surface of these stars are thought to originate in stellar dynamos driven 
in their convective envelopes. There, plasma motions couple with rotation to generate 
global-scale magnetic fields. Though corre lations between the ro tation rate of stars and 
their magnetic activity are observed (e.g.. iPizzolato et al.l l2003) it is at present unclear 
how the stellar dynamo process depends in detail on rotation. 

Motivated by this rich observational landscape, we have explored the effects of more 
rapid rotation on 3-D convection and dynamo action in simulations of stellar convection 
zones. These simulations have been conducted using the anelastic spherical harmonic 
(ASH) code, a tool developed by a team of postdocs and graduate students working with 
Juri Toomre to study global-scal e magnetohydrodynamic convection and dynamo action 
in stellar convection zones (e.g., IClune et all 1 199a iMiesch et al1l200fj iBrun et alJliool 
and contribution by Miesch in these proceedings). 

We began our explorations of convection in rap idly rotating suns by exploring hydro- 
dynamic simulations at a variety of rotation rates (IBrown et al. I l2008l) . These simulations 
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capture the convection zone only, spanning from 0.72i? Q to 0.97i?©, and take solar values 
for luminosity and stratification but the rotation rate is more rapid. The total density 
contrast across such shells is about 25. In those simulations we found that the differential 
rotation generally becomes stronger as the rotation rate increases, while the meridional 
circulations appear to become weaker and multi-celled in both radius and latitude. 

In this paper we review the dynamos we have found in our simulations of more rapidly 
rotating solar-type stars. These wreath-building dynamos form surprisingly organized 
structures in their convection zones ($2j an d some even undergo quasi-cyclic magnetic 
reversals (f3|). Wreath-building dynamos appear throughout the parameter space we have 
surveyed to date (HJ. We close by reflecting on the challenges that lie ahead (fj5]). 



2. A Dynamo with Magnetic Wreaths 

Our first si mulation discussed here, case D3, is of a star rotating three times faster 
than the Sun ( Brown et al.ll2010al) . Vigorous convection in this simulation drives a strong 
differential rotation and achieves sustained dynamo action at relatively low magnetic 
Prandtl number; here Pm=v/i) is 0.5, where v is the viscosity and r\ is the magnetic 
diffusivity. 

The magnetic fields created in this dynamo are organized on global-scales into banded 
wreath-like structures. These are shown for case D3 in Figure QJi. Two such wreaths are 
visible in the equatorial region, spanning the depth of the convection zone and latitudes 
from roughly ±30°. The dominant component of the magnetism is the longitudinal field 
and the two wreaths have opposite polarities. Here the wreath in the northern hemi- 
sphere has negative polarity while the wreath in the southern hemisphere is positive 
in sense. The wreaths are not isolated flux structures; instead, magnetic fields meander 
in and out of each wreath, connecting them across the equator and to higher latitudes 
(Fig. [TJ)) . The lack of visible magnetism in the polar regions reflects the relatively low 
magnetic Reynolds number associated with the convection (average fluctuating Rm' is 
roughly 50 at mid-convection zone). 

It has been a great surprise that such structures can exist in the convection zone of this 
simulation. Generally, it has been expected that convection should shred such structures 
or pump them downwards into a stable tachocline at the base of the convection zone. 
Here the entire domain is convectively unstable, and no such tachocline is present. The 
wreaths persist for long intervals in time, with the mean (longitudinally averaged) mag- 
netic fields showing relatively little variation in time. The convection leaves its imprint 
on the wreaths, with the strongest downflows dragging the field towards the bottom of 
the convection zone. This is visible in the distinct waviness apparent in Figure [TJ On the 
poleward edges the wreaths are wound up into the vortical convection there, and this 
appears to play an important role in regenerating the poloidal field. 

The magnetic wreaths are built by both the global-scale differential rotation and by 
the emf arising from correlations in the turbulent convection. Generally, the mean longi- 
tudinal magnetic field (B^) in the wreaths is generated by the fi-effect: the stretching of 
mean poloidal field by the shear of differential rotation into mean toroidal field. Produc- 
tion of {Btf,) by the differential rotation is balanced by turbulent shear and advection, 
and by ohmic diffusion on the largest scales. 

The mean poloidal field is generated by the turbulent emf 3pi = (u' x £?'), where the 
fluctuating velocity is u' = u — (u) and the fluctuating magnetic fields are B' — B — (B). 
In these simulations, E-pi is generally strongest at the poleward edge of the wreaths, 
centered at approximately ±20° latitude, whereas the 51-effect and (B^) peak at roughly 
±15° latitude. This spatial offset between E-p\ and (B$) means that the turbulent emf is 
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Figure 1. Magnetic wreaths in case D3. (a) Full volume rendering of magnetic wreaths, show- 
ing entire simulation. Lines trace the vector magnetic field with color denoting amplitude and 
polarity of the longitudinal magnetic field (red or light tones, positive; blue or darker tones, 
negative). Rather than being simple flux surfaces, magnetic fields thread in and out of each 
wreath, connecting the wreaths across the equator and linking them to the polar regions, (b) 
Zoom in view of region A showing cross-equatorial connectivity. 

not generally well represented by a simple a-effect description, e.g., 

E F1 = (u' xB')\^^a(B^) (2.1) 

when a is a scalar quantity. This is true even when a is estimated from the kinetic and 
magnetic helicities present in the simulation. More sophisticated mean-field models may 
do much better at matching the observed emf E-pi, and other terms in the mean- field 
expansion may play a significant role; in particular, the gradient of (B^) is large on the 
poleward edges of the wreaths where E-pi is significant. 



3. A Cyclic Dynamo in a Stellar Convection Zone 

We turn now to a more rapidly rotating dy namo simulation, case D5, rotating five 
times faster than the Sun (jBrown et al]|2010bl ). As in case D3, strong global-scale mag- 
netic wreaths are built in the convection zone. Now however, the wreaths begin to show 
significant time-variation and undergo quasi-regular polarity reversals. 

One such reversal is illustrated in Figure [3J Before the reversal (Fig. rjji) , the wreaths 
look similar to those found in case D3, though here magnetism permeates the entire 
convection zone, including the polar regions where relic wreaths from the previous re- 
versal are visible. The equatorial region shows significantly more connectivity and large 
fluctuations of B^, with small knots of alternating polarity visible throughout. This cross- 
equatorial connectivity appears to play an important role in the reversal process. The 
magnetic fields built in case D5 attain somewhat larger amplitudes than those realized 
in case D3: here at mid-convection zone B^, can reach ±40 kG, while in case D3 the peak 
amplitudes were closer to ±26 kG. 

During the reversal (Fig. [2^), new wreaths of opposite polarity form near the equator 
while the old wreaths propagate towards the poles. This poleward propagation appears 
to be a combination of a nonlinear dynamo wave, arising from systematic spatial offsets 
between the generation terms for mean poloidal and toroidal magnetic field, and possibly 
a poleward-slip instability arising from magnetic stresses within the wreaths. After a re- 
versal (Fig. ^jp) the new magnetic wreaths grow in strength and dominate the equatorial 
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Figure 2. Global-scale magnetic reversal in case D5. (a) Magnetic wreaths shown in field 
line tracing shortly before a polarity reversal, with positive polarity wreath above equator and 
negative polarity below. Volume shown spans slightly more than a full hemisphere with both 
polar caps visible where relic wreaths from the previous cycle remain visible. (6) During a 
reversal, new wreaths with opposite polarity form at the equator while old wreaths migrate 
toward the poles, (c) When the reversal completes, the polarity of the wreaths have flipped, 
with negative polarity wreath above the equator and positive below. Cancellation occurs in the 
polar regions and the old wreaths slowly dissipate. Times of snapshots are labeled. 

region. In the polar regions the wreaths from the previous cycle begin to dissipate, recon- 
necting with the pre-existing flux there and being shredded by the turbulent convection 
there. 

Generally, the dynamo processes identified in case D3 serve to build the wreaths of 
case D5. As in case D3, before a reversal the turbulent emf Efi that contributes to the 
mean poloidal field is generally largest on the poleward edge of the wreaths at latitudes 
above ±20°, while the fi-effect and (B^) peak at roughly latitude ±15°. During a reversal 
both E-pi and the production of (B$) associated with the fi-effect surf on the poleward 
edge of the wreaths as those structures move poleward. This systematic phase shift 
appears to contribute to that propagation. 

4. Wreath-building Dynamos 

The two simulations we have explored here, cases D3 and D5, are part of a much larger 
family of simulations that we have conducted exploring convection and dynamo action in 
younger suns. The properties of this broad family are summarized in FigureOi. Indicated 
here are 26 simulations at rotation rates ranging from 0.5 Qq to 15 Qq. At individual 
rotation rates (e.g., 3Q@), further simulations explore the effects of lower magnetic diffu- 
sivity 77 and hence higher magnetic Reynolds numbers. Some of these follow a path where 
the magnetic Prandtl number Pm is fixed at 0.5 (triangles) while others sample up to 
Pm=4 (diamonds). The most turbulent simulations have fluctuating magnetic Reynolds 
numbers of about 500 at mid-convection zone. Wreath-building dynamos are achieved 
in most simulations (17), though a smaller number do not successfully regenerate their 
mean poloidal fields (9, indicated with crosses). Very approximate regimes of dynamo 
behavior are indicated, based on the time variations shown by the different classes of 
dynamos. 

Near the onset of wreath-building dynamo action we generally find little time variation 
in the axisymmetric magnetic fields associated with the wreaths. This is illustrated for 
case D3 in Figure [3f>, showing the mean longitudinal (B^) at mid-convection zone over 
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Figure 3. Wreath-building dynamos, (a) Parameter space showing variety of wreath-build- 
ing dynamos currently explored. Magnetic diffusivity rj and rotation rate SI are shown for dy- 
namo simulations at rotation rates ranging sampling O.5-15f20, with very approximate dynamo 
regimes shown. In some regions, magnetic Reynolds numbers are too low to sustain dynamo 
action, while in other regions persistent magnetic wreaths form which do not show evidence 
for cycles. At higher magnetic Reynolds numbers (occurring here at low 77 or high SI), wreaths 
typically undergo quasi-cyclic reversals. At the highest rotation rates the Lorentz force can 
substantially modify the differential rotation, but dynamo action is still achieved. Cases with 
question marks show significant time-variation but have not been computed for long enough 
to definitively establish cyclic behavior. (6) Time-latitude plot of mean (axisymmetric) B$ at 
mid- convection zon e in persistent case D3 (|Brown et al.ll2010al ). (c) Cyclic case D5 shown for 
same span of time l| Brown et a l. 2010b). Three reversals are visible here, occuring on roughly 
1500 day periods. Times of snapshots shown in Figure [2] are indicated. 

an interval of about 4000 days. Though small variations are visible on a roughly 500 day 
timescale, the two wreaths retain their polarities for the entire time simulated (more than 
20,000 days), which is significantly longer than the convective overturn time (roughly 10- 
30 days), the rotation period (9.3 days), or the ohmic diffusion time (about 1300 days 
at mid-convection zone). We refer to the dynamos in this regime as persistent wreath- 
builders. 

Generally, we find that wreath-building dynamos begin to show large time dependence 
as the magnetic diffusivity 77 decreases and as the rotation rate S\ increases. In many 
cases this leads to quasi-regular global-scale reversals of magnetic polarity, as discussed 
for case D5 in $3] We illustrate three of these cycles occuring in case D5 in Figure [3f>. In 
this simulation, reversals occur with a roughly 1500 day timescale, though during some 
intervals the dynamo can fall into other states. For reference, the ohmic diffusion time in 
this simulation is about 1800 days, while the rotation period is 5.6 days. Similar cycles 
occur in other simulations, but the period of reversals varies with both 77 and Q. Cycles 
appear to become shorter as r\ decreases, opposite to what might be expected if the ohmic 
time determined the cycle period. The dependence of cycle period on SI is less certain. At 
present, determining why cycles are realized in many of these dynamos remains difficult. 
The phenomena appears to be at least partially linked to the magnetic Reynolds number 
of the differential rotation and possibly to that of the fluctuating convection. 

In these wreath-building dynamos the major reservoir of kinetic energy that feeds the 
generation of magnetism is the axisymmetric differential rotation, and this global-scale 
shear is strongly reduced in the dynamo simulations. Individual convective structures 
are largely unaffected by the magnetic wreaths except when the fields reach very large 
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Figure 4. Nests of convection in dynamo case D10L. Convective patterns are shown by snap- 
shots in Mollweide projection of the radial velocities v r (a) near the surface (0.95 Rq) and (b) at 
mid- convection zone (0.85 Rq). Two active nests of convection are clearly visible, (c) Near the 
surface, the radial magnetic field B r is concentrated in the stronger nest. This patch of radial 
field propagates with the nest, (d) At mid-convection zone the active nests leave their imprint 
on the magnetic wreaths, shown here by longitudinal field B^. 



amplitudes; in case D5 this occurs when exceeds values of roughly 35 kG at mid- 
convection zone. At the highest rotation rates the Lorentz force of the axisymmctric 
magnetic fields becomes strong enough to substantially modify the differential rotation, 
largely wiping out the latitudinal and radial shear (e.g., cases D10 and D15 in Fig. [3]i). 
In these cases, wreath-like structures can still form though they typically have more 
complex structure and are less axisymmetric. 

Some of these rapidly rotating dynamos have held further surprises. One such simula- 
tion, case D10L rotating ten times faster than the Sun, is shown in Figure |H Magnetic 
wreaths form in this dynamo, but so do strong localized nests of convection. Two such 
nests are visible near the equatorial region in the convective patterns of radial velocity 
(Fig.|3h, b). At certain longitudes the convection is significantly stronger, and these nests 
span the convection zone. The nests leave their imprint on the magnetism, concentrat- 
ing the radial magnetic field B r into patches near the surface (Fig. Hh) and interacting 
strongly with the toroidal field B^ at depth (Fig. 0]i) . These active nests appear to be dis- 
tinct dynamical structures, and they propagate at a different rate than either individual 
convective cells or the local differential rotation. If such nests persist in stellar convec- 
tion zones, they are likely to strongly affect inferences of the global-scale m agnetic field 



We h ave found active nests in rapidly rotating hydrodynamic simulations ([Brown et al 



2008), but to our knowledge this is the first time such structures have been found in a 
fully saturated dynamo state. 

We have found wreath-like structures in simulations rotating at the solar rotation rate 
as well (one example appears in the contribution by Miesch, these proceedings). Achieving 
such structures in solar simulations is challenging, as the angular velocity contrast in 
the Sun is smaller than that realized in the rapidly rotating dynamos. Generally the 
solar simulations require low values of r\ to build wreaths, which in turn calls for high 
resolutions and that exacts a large computational cost. 



Wreath-building dynamos 



7 



The magnetic boundary conditions adopted near the base of the convection zone also 
play an important role: here we have explored simulations with perfectl y conducting bot- 
tom boundaries. Previous simulations of the solar dynamo (case M3. iBrun et al.l 120041) 
used a potential field bottom boundary and did not find magnetic wreaths like those 
shown here. When the boundary conditions are changed in case M3 to the perfect con- 
ducting boundaries used here, the axisymmetric magnetic fields grow in strength, become 
more organized and form wreaths near the base of the convection zone. A variant on that 
solar dynamo simulation is shown in the contribution by Miesch in these proceedings. 
Conversely, magnetic wreaths are diffic ult to achieve in the rapidly rotating dynamos 
when we use the boundary conditions o f B run et al. (2004). We take some comfort in the 
fact that wreaths continue to form in simulations which include a portion of the stable 
radiative zone within the domain. In the simulations with model tachoclines, the wreaths 
of magnetism are essentially unmodified from those found in the simulations which only 
capture the convection zone. This suggests that a perfectly conducting bottom boundary 
may mimic the presence of a highly-conductive radiative zone below the convection zone 
better than a potential field extrapolation does. 



5. Overview 

Stellar convection spans a vast range of spatial and temporal scales which remain 
well beyond the grasp of direct numerical simulation, and we remain humbled by the 
complexities posed by highly turbulent convection on global-scales in rotating, stratified, 
and magnetized plasmas. Stellar dynamo studies must drastically simplify the physics of 
the stellar interior: as an example, molecular values of rj in the solar convection zone range 
from roughly 10 2 -10 5 cm 2 /s as one moves from the tachocline to the near surface regions, 
while the molecular viscosity v is of order 10 cm 2 /s there. In contrast, our simulations 
employ values of r\ and v that are of order 10 12 cm 2 /s; this large value is more similar to 
simple estimates of turbulent diffusion associated with granulation at the surface where 
v t ~ VtLt ~ 10 11 cm 2 /s given Vt ~ 1 km/s and L% ~ 1 Mm. Despite this daunting 
separation in parameter space, it is striking that coherent magnetic structures can arise 
at all in the midst of turbulent convection. We find the combination of global-scale spatial 
organization and cyclic behavior fascinating, as these appear to be the first self-consistent 
3D convective stellar dynamos to achieve such behavior in the bulk of the convection zone. 

A variety of wreath-building dynamos have been found in these simulations of rapidly 
rotating suns: some build persistent wreaths, while others undergo significant time vari- 
ations including quasi-cyclic reversals of global-scale magnetic polarity. The parameter 
space is complex, and some simulations show cycles in one hemisphere but not the other. 
Cyclic cases like case D5 can even wander into and then back out of non-cycling states. 
The role of tachoclines in stellar dynamos remains a matter of great debate. These sim- 
ulations suggest that, at least in rapidly rotating stars, tachoclines may not play as 
crucial a role in the organization and storage of the global-scale magnetic field as in the 
solar dynamo. A major step forward will be to explore simulations that couple wreath- 
building dynamos in the convection zone, through a tachocline of shear, to the stable 
radiative interior below; those simulations are ongoing now. We are also exploring how 
convection and dynamo action may be different in other solar-type stars, including K- 
type dwarfs and the fully convective M-type stars, which present their own mysteries 
(see Browning! 120081 and contribution by Browning, these proceedings). The future of 



3D dynamo simulations is very bright, and cycl i c solutions are beg i nning to appear i n 
a variety of situations (e.g., Ghizaru et al.|[201ot iKapyla, et al. 20101 : Mitra et al.l 20ldh 



These are exciting times indeed for stellar dynamo theorists! 
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Discussion 

T. Rogers: Could the poleward propagation of held be likened to the dynamo wave, 
in which the sign of propagation depended on helicity associated with convection versus 
radiation zones? 



B. Brown: Maybe. Here we clearly see signs that the poleward propagation is partly due 
to Maxwell stresses and a resulting poleward slip, and this phenomena occurs essentially 
unmodified in simulations that include a tachocline and convective penetration into a 
stable radiative zone. But there may be dynamo wave aspects to the reversals as well. 
[Note: subsequent work indicat es that a non- l inear d ynamo wave does play an important 
role in the reversal process; see iBrown et al. (2010b).] 



E. Zweibel: Why is the "failed dynamo" region on the r)-VL plane localized near the 
solar rotation rate (ft — 51©)? (see Figure ISJi) 

B. Brown: On either side of the solar rotation rate, the latitudinal shear is strong. 
Simulations spinning slower than the Sun have strong anti-solar differential rotation 
(retrograde equators, prograde poles) while the more rapidly rotating simulations have 
stronger solar- like differential rotation (prograde equators, retrograde poles). Both cases 
lead to higher magnetic Reynolds numbers associated with the latitudinal shear. Simu- 
lations at the solar rotation rate have weaker differential rotation, and need lower values 
of r\ to achieve dynamo action in these wreath-building dynamos. 



